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It is shown that the bond-decorated Ising model is a realistic model for certain real magnetic 
compounds containing lanthanide ions. The lanthanide ion plays the role of Ising spin. The re- 
quired conditions on the crystal-field spectrum of the lanthanide ion for the model to be valid are 
discussed and found to be in agreement with several recent ab initio calculations on Dy 3+ centers. 
Similarities and differences between the spectra of the simple Ising chain and the decorated Ising 
chain are discussed and illustrated, with attention to level crossings in a magnetic field. The mag- 
netic properties of two actual examples (a [DyCuMoCuJoo chain and a Dy 4 Cr 4 ring) are obtained 
by a transfer-matrix solution of the decorated Ising model, g-factors of the metal ions are directly 
imported from ab initio results, while exchange coupling constants are fitted to experiment. Agree- 
ment with experiment is found to be satisfactory, provided one includes a correction (from ab initio 
results) for susceptibility and magnetization to account for the presence of excited Kramers doublets 
on Dy 3+ . 
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FIG. 1. The Ising-Heisenberg chain discussed in Ref. |3] (show- 
ing only two unit cells). The bonds of the Ising chain (1-1'- 
1"- . . . ) are decorated with Heisenberg dimers ((2-3), (2'-3'), 
. . . ) . The partition function of this chain is exactly solvable. 



I. INTRODUCTION 

In a 1959 paper M. E. Fisher introduced the general 
bond-decorated Ising model as one example of a set of ex- 
actly solvable transformations of spin-1/2 Ising models. 1 
A bond-decorated Ising model has an "arbitrary statis- 
tical mechanical system" inserted in every original Ising 
bond. The partition function of this decorated model is 
related to the partition function of the original or bare 
Ising model by the addition of a prefactor and a renor- 
malization of the coupling constants and magnetic mo- 
ments (the Ising model is supposed to be in a parallel 
magnetic field). 2 Knowledge of the partition function of 
a given Ising model thus allows one to obtain the parti- 
tion function of any bond-decorated version of that Ising 
model. 

More recently, Strecka and Jascur have used the 
method of bond-decoration to investigate the thermo- 
dynamics of mixed Ising-Heisenberg chains in paral- 
lel magnetic fields where the decorating unit is a spin 
dimer or trimer with anisotropic Heisenberg coupling (see 
Fig.JlJP^The partition function of these chains is readily 
obtained from the known partition function of the Ising 
chain and the energy levels of the decorating unit. There- 
fore, they could calculate exact magnetic properties and 
theoretically show, for example, the existence of magne- 
tization plateaus in certain bond- alternating chains. 



The convenience of the decorated Ising chain as a the- 
oretical model for spin chains derives from the relative 
ease with which exact solutions are obtained, in contrast 
for example with the pure Heisenberg chain, for which 
no exact partition function has been found. Up to now, 
this property of solvability has been the prime motive 
for the study of these chains in the literature. Indeed, in 
Ref. [3] the decorated Ising chain was considered as a sub- 
stitute for the intractable Heisenberg model, and in Ref. 4 
the principle reason for introducing Ising bonds — to re- 
place the more reasonable Heisenberg bonds — in a chain 
of Cu 2+ ions was the desire to obtain a solvable model. 
This approach can be applied to any type of Heisenberg 
chain with a repeating unit: replace enough — preferably 
ferromagnetic — Heisenberg bonds by Ising bonds to ob- 
tain a decorated Ising chain that is solved easily and, 
in some cases, exhibits thermodynamic properties that 
are q uali tatively comparable with those of the original 
chain.EH 

However, the role of the decorated Ising chain in the 
field of one-dimensional magnetism is not confined to 
that of a simplified model of realistic quantum spin 
chains. In this paper we show that some new molecu- 
lar rings and chains are real examples of decorated Ising 
systems. Concretely, we treat a [DyCuMoCu] infinite 
chairPand a (DyCr)4 tetrameric ringP These compounds 
were recently synthesized in the course of the ongoing 
synthetic efforts to make new and better single-chain 
magnets (SCMs) and single- molecule magnets (SMMs), 
whose characteristic property is a blocking or slow relax- 
ation of magnetization at low temperatures. 9 We will not 
be concerned here with these dynamical aspects of their 
magnetism, but only with their static magnetic prop- 
erties. A necessary property of SMMs and SCMs is a 
magnetic anisotropy. One line of approach is to intro- 
duce anisotropy by means of lanthanide ions, wheth er or 
not in combination with transition metal ionsP^El The 
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\ I 1 2. Part of the decorated Ising chain. In each Ising 
bond is inserted an arbitrary statistical mechanical systerrP, 
called decorating unit, that interacts only with the two Ising 
spins Si at the vertices of the bond. The Ising spin variables 
commute with the Hamiltonian and have definite values in 
the eigenstates of the chain. 



two compounds considered here are products of this ap- 
proach, with dysprosium as lanthanide ion. 

The Dy 3+ ion plays a crucial role in these systems; 
the nature of the ground state of this ion in its ligand 
environment determines whether the system is a deco- 
rated Ising system or not and consequently, whether its 
partition function is exactly solvable or not. The ground 
Kramers doublet of Dy 3+ must have complete uniaxial 
magnetic anisotropy^ and must be separated from ex- 
cited Kramers doublets by an amount that is large com- 
pared with the exchange coupling (typically, this sepa- 
ration must be 100 cm" 1 or more) P^l The required in- 
formation on the ground and excited doublets of the 
Dy 3+ monomer can be derived from ab initio calculations 
on the mo nom er complex, isolated from the polynuclear 
compound.^ 

The [DyCuMoCu] chain and the (DyCr)4 ring are 
shown to be decorated Ising chains in an arbitrarily di- 
rected magnetic field. The magnetic properties, in partic- 
ular powder magnetization and susceptibility, are calcu- 
lated with the help of the transfer-matrix method, which 
is a bit more general and convenient for numerical compu- 
tation than the renormalization of the Ising parameters, 
which was used by Fisher and by Strecka and Jascur. 
The results compare well with experiment and allow to 
determine values for the exchange coupling constants. 

The excited crystal field Kramers doublets (or Stark 
levels) of Dy 3+ are not included in the decorated Ising 
model. Because of their relatively low energetic position, 
these Kramers doublets can have a non-negligible contri- 
bution to the magnetic properties of the chain. We ac- 
count for this in a first approximation by adding this con- 
tribution, as calculated ab initio for the monomeric Dy 3+ 
complex, to the results of the decorated Ising model. 



II. THEORETICAL BACKGROUND 

A. Decorated Ising chain and transfer-matrix 
solution 

The decorated Ising chain may be divided into units 
delimited by the Ising spins, as in Fig. [2j The Hamil- 
tonian of the decorated Ising chain (or ring) of length n 
may accordingly be written as follows: 

H = h 1 (s 1 , s 2 )+h 2 (s 2 , s 3 )+h 3 (s 3 , s 4 ) + . . .+h n (s n , s n+1 ), 



where the subhamiltonians hi correspond to units of the 
chain. If the chain is closed into a ring, periodic bound- 
ary conditions apply by identifying the last with the first 
spin: s n+ i = s\. The Ising spins Sj commute with the 
subhamiltonians and the subhamiltonians commute with 
each other: 



[Sj, hj] =0 and 



[hi,hj] = 0. 



(2) 



It is further assumed that there is no direct interaction 
between the decorating units themselves, i.e., two differ- 
ent hi have no variables in common, except for an Si when 
they are neighbors. In writing Eq. (JlJ we made use of 
the fact that the Ising spins are conserved variables and 
may be considered as parameters rather than operators. 
In the case of spin- 1/2, they take on the values +1/2 and 
— 1/2, so that to each decorating unit there correspond 
four different subhamiltonians hi(si, Sj+i)PI 

Since there is no direct interaction between the dec- 
orating units of the chain, the subhamiltonians hi in 
Eq. ([I]) work — for a given set of Sj values — on disjunct 
spaces. An eigenfunction of H is then simply a direct 
product of n independent eigenfunctions, one of each hi- 
The corresponding energy is: 



E = e lkl {s 1 ,s 2 ) +£2k 2 (s 2 ,s 3 ) + e 3k3 (s 3l s 4 ) + 

+ £nfe„ Sn+l), 



(3) 



where Sik t (s*, Sj+i) is the fcj-th energy level of hi{si, Sj+i). 

Following Ref. [T]we write the partition function of one 
decorating unit for fixed s$ and Sj+i on the bond vertices: 

ipi(si,s l+1 ) = exp[-/3£ik(sj, s i+ i)], (4) 



where (3 = 1/fcT. The total partition function for the 
chain is then given by 

This is the well-known form of a transfer-matrix 
solution!^ Each decorated bond (or pair of neighboring 
Ising spins) is represented by a transfer matrix 2$ whose 
elements are the ipi( s i> s i+i)> the values of s, labeling the 
rows and the values of Sj+i labeling the columns. Explic- 
itly for the spin-1/2 decorated Ising chain: 



i_ UH>+I) M-l-l) 



(6) 



Since most applications deal either with rings or very 
long chains, the boundaries can be identified (s„+i = Si) 
and expression (|5| can be written as the trace of a matrix 
product: 



Z n = Tr(T 1 T 2 ...T n ). 



(7) 



(1) 



This is the most general expression for the partition func- 
tion of a decorated Ising chain. If the chain has trans- 
lational symmetry the partition function is expressed in 
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terms of eigenvalues.^ Suppose that the chain is a repe- 
tition of identical decorated Ising bonds, then the Tj are 
all the same, so that 



Z„ = Tr(T 1 ") = ^A^ 



(8) 



where the A^ are the eigenvalues of T\. According 
to Perron's theorem on positive matrices, the eigen- 
value Ao with largest modulus is real, positive and 
nondegenerateP^l This yields the simple but exact result 
for the free energy per unit cell of the infinite chain: 



In Z„ 



f = —kT lim = — fcTlnAo- 



(9) 



If a unit cell of the chain is spanned by p decorated 
bonds instead of one then it is only necessary to com- 
bine p transfer matrices into one new transfer matrix 
T = T{F 2 . . . TL, with largest eigenvalue A , to be used in 
Eqs. ^ and ([9]), where n is to be replaced by n/p. This 
situation arises, for example, when the local easy axes on 
the Ising ions are not parallel but canted with respect to 
each other. 

The solution in terms of transfer matrices, Eq. Q, 
is not limited to decorated spin-1/2 Ising chains. It is 
valid for chains having Ising spins of any multiplicity or 
a combination of Ising spins with different multiplicities, 
in which case the dimension of the transfer matrix is dif- 
ferent from two by two. Another advantage of the trans- 
fer matrix method is that it can be readily extended to 
include ncxt-ncarcst-neighbor bonds between the Ising 
spins. To this end, the transfer matrix has to be en- 
larged so that it does not jump from one Ising spin to 
the next, but from one pair of Ising spins to the next 
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pair. 



B. Dy 3+ as Ising spin 

We are interested in this paper in molecular chains or 
rings that can be described by a decorated Ising model. 
This means that their low-energy spectrum can be mod- 
eled to satisfactory accuracy by an effective Hamiltonian 



that has the properties described in Section II A It has 
been noted there that the composition of the decorating 
unit is basically arbitrary and we need therefore not con- 
sider the properties of that part. Instead, our attention 
goes here to the molecular realization of the Ising spin [sj 
in Eq. Q]. We focus on Dy 3+ because this ion is used 
in the two examples studied below, but the discussion 
applies equally well to several other trivalent lanthanides 
and also to some transition metal ions. 

It is known that lanthanide ions in a coordination envi- 
ronment are often well described by crystal field theory 
applied to the ground \LSJ) level. One assumes that 
J, L, and S remain good quantum numbers. Dy 3+ is a 
Kramers ion that belongs to the second half of the lan- 
thanide series, whose ground level is 6 H 15 y 2 (/ 9 ), with 



associated Lande factor g = 4/3. This multiplet splits 
into eight Kramers doublets by the crystal field perturba- 
tion (except for high-symmetric environments belonging 
to the T, O, or K point groups, which split the multiplet 
in less than eight levels). It will be useful to view each 
Kramers doublet as an effective spin-1/2 with its own g- 
factors (3 in number) and corresponding magnetic axes. 
For example, take the Kramers doublet \Mj = ±15/2) 
(of the 6 H 15 ^ 2 level), quantized with respect to the z axis. 
Its g factors are g z — 2(4/3)(15/2) = 20 and g x = g y = 0. 

If the action of the crystal field on Dy 3+ is such that 
the lowest Kramers doublet is separated from the next 
one by an energy that is large compared to the energy 
of interaction with the magnetic field and the exchange 
interaction with neighboring ions, we can omit all ex- 
cited Kramers doublets from the Hamiltonian and keep 
only the lowest doublet. In this way Dy 3+ is described 
by a spin of 1/2 and every interaction in which it takes 
part enters the Hamiltonian as a linear combination of 
the three spin operators s x , s y , and s^P^If we now want 
Dy 3+ to be an Ising spin, defined by the first commuta- 
tion relation in Eq. ([2J), it is clear that only one of these 
spin components, say s z or simply s, may actually appear 
in the Hamiltonian. In other words, there must be no in- 
teraction that creates an off-diagonal matrix element be- 
tween the two components of the Kramers doublet. This 
can be shown to be true with high accuracy if the lowest 
Kramers doublet is \Mj = ±15/2). 

The two interactions of importance here are the Zee- 
man interaction with the magnetic field and the ex- 
change interaction with other magnetic ions. The Zee- 
man Hamiltonian follows directly from the g factors of 
the Kramers doublet (vide supra), and is 



B z 



(10) 



where B z is the z component of the applied magnetic 
field. Note that the field may be applied in any direction, 
but it is only the z component that interacts with the 
Kramers doublet because g x = g y = 0. The vanishing 
of g x and g y in \Mj = ±15/2) follows from the selection 
rule stating that a vector operator (the magnetic moment 
in this case) cannot connect states for which Mj differs 
by more than one unit. We will sometimes refer to the 
z axis as the anisotropy axis, to stress that it is the only 
magnetic axis with nonvanishing g factor. 

To evaluate the effect of exchange interaction, we 
must first take a closer look at the composition of the 
Kramers doublet. In terms of the Russel-Saunders states 
\M S )\M L ) we have 



15/2) 
15/2) 



■5/2)| + 5) 
5/2)|-5). 



(11) 



In a basic (super)exchange process between two mag- 
netic centers, one electron of each center takes part. If 
we look at one center, the process removes an electron 
with certain spin projection (up or down) and puts it 
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back on the center either with the same or with oppo- 
site spin projection.^ This gives rise to the selection rule 
AMs = 0, ±1. If the exchange interaction is to connect 
both components of the Kramers doublet in Eq. (11), at 



necessarily of the form 



least five successive processes are needed, for AMs = ±5- 
In other words, not one but five electron spins have to be 
flipped to connect | + 15/2) with | — 15/2). If the basic 
exchange process (i.e. the one for which AMs = 0, ±1) 
occurs in, say, kth ordei^ of perturbation theory then 
an off-diagonal matrix element between the two compo- 
nents can only appear in (5fc)th order of perturbation 
theory. It is therefore reasonable to assume that the off- 
diagonal matrix element is negligibly small compared to 
the diagonal matrix elements (AMs — 0) so that the ef- 
fect of exchange interaction on the Kramers doublet is 
accurately described by the s z spin operator only. 

Note that we derived the selection rule for exchange 
interaction on the basis of the spin quantum number Ms 
only, without paying attention to the angular momentum 
quantum number Ml, although Ml changes even more 
than Ms between the states of the Kramers doublet (111. 



The existence and precise form of a selection rule for 
Ml depends on the spatial symmetry of the exchange 
problem under consideration. AMf, is therefore not as 
useful as AMs for predicting the vanishing of certain 
matrix elements of exchange interaction. Note however 
that, even in the lowest symmetry, there is a maximum 
to the amount that Ml can change in the basic exchange 
process described in the previous paragraph: within the 
/ orbitals, a one-electron process can bring about at most 
a change of AMl = ±6. So we would have, in general, 
that a basic exchange process can induce the following 
changes in a lanthanide state: 



AM 5 = 0, ±1 
|AMl| < 6. 



(12) 



Thus at least two steps of this kind are needed to bridge 
AMl = ±10, but at least five are needed to bridge 
AMs = ±5. So in this case, the selection rule of Ms 
gives the stronger result and leads to the conclusions — 
reached in the last paragraph — on the matrix elements 



of exchange interaction in the Kramers doublet ( 11 ) 



We can now derive the precise form of that part of the 
effective Hamiltonian that refers to the exchange inter- 
action between a Dy 3+ ion [with ground state (111] and 
another magnetic center. We consider two cases, which 



we shall encounter in the examples in Sections III and 



in the first case the other center is another Dy 



3+ 



IV 

ion; in the second case the other center is an ion with 
an isotropic spin moment S. Consider first the exchange 
interaction with another Dy 3+ ion. We assume that the 
second Dy 3+ has the same property of uniaxiality as the 
first one and that it also shares all other relevant prop- 
erties discussed in the previous paragraphs. Then both 
ions are represented by a spin-1/2 doublet with a local 
anisotropy axis z, (i = 1 , 2) and we already know that 
the effect of exchange interaction in each doublet is pro- 
portional to s Zi . Therefore the exchange Hamiltonian is 



(13) 



where it is understood that the first spin belongs to ion 
1 and the second to ion 2. Note that Zi and Z2 need not 
be parallel with each other. 

As a second case, consider the interaction between 
Dy 3+ (anisotropy axis Zi) and an isotropic spin S2. 
The latter may typically be a transition metal ion with 
quenched orbital momentum. We found that the ex- 
change processes that contribute do not change the spin 
projection on the Dy 3+ ion (AMg 1 = 0, quantization 
axis Zi), and that this result is independent of the ex- 
change partner. It is also known that every exchange 
process commutes with the total spin (the matrix el- 
ements involved are spin-independent matrix elements 
of kinetic, potential, and Coulomb energy^^), so that 
AM Sl + AM Sa = 0. It follows then that AMs 2 = 0, or, 
the exchange Hamiltonian commutes with the Z\ compo- 
nent of S2- The simplest expression compatible with this 
requirement is 



J s Zl S Zl , 



(14) 



where s naturally represents Dy and S represents S2. 
The interaction is of Ising form with the anisotropy axis 
of Dy 3+ as Ising axis. Note that, when S2 > 1/2, higher 
powers of S Zl may enter the Hamiltonian. Consider- 
ing exchange interaction as a perturbation however, one 
can usually assume that the lowest-order contribution, 



Eq. (14), is the leading term. 



There are other cases conceivable, for example dipole- 
dipole interaction between the moment of Dy 3+ and a 
neighboring moment. One will always find, as above, that 
the Hamiltonian is a product of s z (belonging to Dy 3+ ) 
and a part that belongs to the other ion and whose form 
depends on the kind of the other ion and on the details 
of the interaction. 

We have now obtained that a Dy 3+ ion, if its ground 
state is \Mj = ±15/2), fairly well separated from excited 
states, interacts with the magnetic field and with neigh- 
boring ions as an Ising spin-1/2, in the sense that the in- 
teraction is always proportional to s z , as expressed in the 
Eqs. (10), (13), and (14). This means that a chain-like 
molecular structure having Dy + ions of this kind at reg- 
ular positions in the chain would meet the requirements 
of a decorated Ising chain, given in part by Eqs. ([I]) and 
([2]). It remains of course to be shown that |Mj = ±15/2) 
can indeed be the ground state of a coordinated Dy 3+ 
ion in a polynuclear complex. 

At first sight, this seems rather unlikely. \Mj = 
±15/2) is an eigenstate of cylindrical symmetry. Within 
lanthanide f n states, the crystal field is effectively 
cylindrical if there is, at least, an eightfold rotation 
axis (Cs) or rotation-inversion axis (Sg)!^ 5*8 sym- 
metry has been obtained, for example, in mononu- 
clear bis(phtalocyaninato) sandwich complexes of the 
lanthanides. 21 Even when such high symmetry is at- 
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tained, the ground state is not necessarily the cylindri- 
cal doublet with highest \Mj\ value Apart from that, 
the symmetry of the coordination sphere of a lanthanidc 
ion in a polynuclear, possibly heterometallic, complex 
or chain is usually much lower or even completely ab- 
sent. Such is the case for the two examples considered 
in this paper. On the basis of symmetry alone, there is 
thus no reason to expect \Mj = ±15/2) to be an eigen- 
state, let alone the ground state. Nevertheless, recent 
ab initio calculations have revealed the unexpected re- 
sult that the ground state of several low-symme try com - 
plexes of Dy 3+ is very close to \Mj = ±15/2) 
They used the multiconfigurational, wavefunction-based 
CASSCF/RASSI-SO method, to obtain accurate wave- 
functions for several of the lowest Kramers doublets. Cal- 
culation of the principal g factors of these states gives an 
indication of their composition. It was found in several 
cases that the ground doublet has g z close to, but lower 
than 20, and g x and g y close to 0. This corresponds to a 
doublet mainly composed of \Mj — ±15/2). 

With evidence of ab initio calculations it is thus pos- 
sible to identify Dy 3+ in certain coordination environ- 
ments as an Ising spin-1/2 (to a good approximation). 
We will use this information to identify the compounds 
in Sections . 



To cone 



III and IV as decorated Ising chains (or rings) 



ude this section we remark that one can- 
not deduce from the vanishing of two g factors alone 
that a Kramers doublet will behave as an Ising spin. 
It will, of course, in its interaction with the magnetic 
field [Eq. (uOj], but this will not, in general, be true 
for the exchange interaction. Take again the example 
of Dy 3+ , supposing the ground state is the Kramers 
doublet \Mj = ±7/2). It is perfectly uniaxial because 
9x = 9y — and g z = 9|. A closer look at the expan- 
sion of \Mj = ±7/2) in terms of the Russel-Saunders 
states \Ml)\Ms) shows, however, that the selection rules 
in Eq. \12\ permit a matrix element to exist between 
1+7/2) and | — 7/2), therefore introducing non-Ising terms 
(i.e., s x and s y ) in the same order of perturbation theory 
as the Ising term in the exchange Hamiltonian. 



C. Nature of the eigenstates and level crossings 

There exists a similarity between the spectra of the 
decorated and undecorated Ising chains at which wc 
want to take a closer look here. We suppose infinite, 
periodic chains, or periodic, even-membered rings (In 
odd-membered rings spin-frustration complicates the pic- 
ture). We also limit ourselves to chains with Ising spins 
of 1/2. 

The undecorated, or simple, Ising chain in a magnetic 
field B parallel with the z axis is given by Eq. ([I]) and 



hi = -JsiS i+1 - n B gSiB, 



(15) 



where Si is the z component of Sj (si = ±1/2). We as- 
sume, without loss of generality, B > 0. The eigenstates 
of the chain are spin configurations like (ttiti ■■■)■> etc. 



Of the 2" eigenstates only two distinct ones can be the 
ground state: the ferromagnetic (F) and the antiferro- 
magnetic (AF): 

F:(tttt-..) if J>0 or J < and ^gB > \ J\, 
AF : (UH ■■■) if J < and fi B gB < \ J\. 

(16) 

When B = time-reversal symmetry makes every state 
degenerate with the state formed by flipping all the spins. 
This degeneracy is meant to be implied in (16 1, where 
only one of two states is shown in each case. Only the two 
AF states remain degenerate when B ^ 0. When J < 
a ground state level crossing occurs from AF to F when 
B is increased. At the point of crossing (/iBgB = \J\) 
the two AF states are degenerate together with all states 
derived from an AF state by flipping one or more down- 
spins up. However, no other state than AF or F can 
be the ground state at any other value of B. Thus the 
ground state of the Ising chain is either F or AF, and 
they are degenerate, together with an infinite number of 
other states, at the crossing point. 

We now decorate the Ising bonds with identical but ar- 
bitrary units to obtain a periodic decorated Ising chain. 
The spectrum is given by Eq. ([3|. Since the chain is 
periodic, the spectrum of the individual units hi is inde- 
pendent of i and the energies may be written as with 
k ranging over the eigenstates of hi . There are four sets of 

£k- £fc(tt)> Efe(t-J-), £fc(4-t)) and £fc(44)> in an obvious no- 
tation. Notice that the eigenstates of this chain can still 
be classified according to the configuration of the Ising 
spins: (ftitl ■ • • ) etc., which follows from the fact that 
all the Si and H form a commuting set of observables. 

An interesting question is whether the same rules hold 
for the ground state of the decorated Ising chain as did 
for the simple Ising chain. The answer is yes; the ground 
state is either F or AF (referring to the Ising spin con- 
figuration) and a crossing between them is possible, with 
the same number and kind of degenerate states as in the 
simple Ising chain. To show this, we have to consider 
only the lowest eigenstate belonging to each of the 2™ 
possible Ising spin configurations. In these states every 
unit is in its lowest possible state for the given orientation 
of the neighboring Ising spins: £i(s,s') (we assume this 
energy to be nondegenerate) , so that the total energy of 
the chain state is 

E = n tt e 1 (tt)+n w .ei(44,)+n W Ei(|t)+»H4. e i(t4'). ( 17 ) 

where n^f denotes the number of pairs of neighboring 
Ising spins that are both spin up, etc. For example, in 
the F configuration in Eq. ( 16 ), riff = n (periodic bound- 
ary conditions are assumed), while in the AF configura- 
tion, Tiff — 0. The eigenstates we have just described, 



with energy (17), are in an obvious one-to-one correspon- 



dence with the eigenstates of the simple Ising chain. The 
ground state is found by minimizing ( 17) with respect to 



G 



the n ss > , under the restrictions 



n n + n u + n it + n n = n 



(18) 



The first relation states that the total number of Ising 
spins (or, equivalently, unit cells) is n. The second re- 
lation follows from the fact that, in a cyclic spin config- 
uration, every J,f pair must eventually be followed by a 
ti pair, possibly after a number of ft pairs. Another re- 
striction is that whenever both riff and are not zero, 
[and by Eq. ( 18 1 also rtn ] must be at least one. Using 
(18) we can rewrite Eq. <JT7J) as 



E : 



: «n( £ i(tt) 

•ra.U-( s i(44-) - 



ei(n)+£iUt)] 



(19) 



MtD+eidt)], 



where we see that only the average [£i(tl) + £i(lt)]/2 of 
the "antiparallel" energies enters the equation. The last 
term is a constant and can be discarded for the purpose 
of relative energy considerations. 

We can now derive the values of riff and for the 
ground state of the chain, keeping in mind that £i(s, s') 
is a function of the magnetic field B. Suppose then, 
first, that B = 0. Time reversal symmetry asserts that 
£ i(tt) = £ i{\X) an( l £i(tl) = £i(4t)< It is simple to see 
that, depending on the relative ordering of £i(tt) an d 
e i(t4-)i E is minimal in the F configuration (riff = n or 
n u = n; n u = riff = ) when £i(tt) < £i(tl) or in 
the AF configuration (riff = = 0; = nff = n/2) 
when £i (tt) > £i (ti) ( we exclude the possibility of equal- 
ity of both energies from the discussion; £i(tt) = £ i(tl) 
would correspond, in the simple Ising chain ( |15[ ), with 
J = 0). When B ^ 0, time reversal symmetry is not op- 
erative, and we have, in general, four different energies 
£i(s,s'). The equation (19) shows that the configura- 



tion that minimizes E is determined by the sign of the 
two terms in round brackets; if both are positive, then 
nff = = (AF configuration); if at least one of them 
is negative, then either nff = n or = n (F configura- 
tion), depending on whether respectively £i(tt) ° r 
is lower. 

Finally, the magnetic field can induce a transition from 
the AF to an F ground state configuration, say with all 
Ising spins up. This happens when 



ei(tt) = i[ei(tlO+eiUt)] 



(20) 



and £i(tt) < At this point, the ground state 

configurations are all those for which {n^ = 0, riff = 
0, 2, 4, . . . , n}, exactly the same as in the simple Ising 
chain. 

We find thus a complete analogy between the simple 
and the decorated Ising chain as far as the ground state 
Ising spin configuration is concerned. The only possible 



configurations are the fully antiferromagnetically aligned 
and the fully ferromagnetically aligned configurations. 
No "intermediate" configuration can be the ground state. 
The only exception is the crossing point between AF and 
F, where there is a high degeneracy of configurations. 
These conclusions are independent of the nature of the 
decorating unit. 

Although the decorated Ising model predicts that the 
AF and F ground states are both doubly degenerate (in 
zero field), this degeneracy is not a result of the spatial 
symmetry: in the cyclic group C n , the two AF compo- 
nents combine into irreducible representations (irreps) A 
and B, while the two F components transform as two A 
irreps. Introduction of neglected terms in the Hamilto- 
nian, that destroy the Ising property, could split these 
ground state components. 

The decorated Ising chain affords two new kinds of 
ground state level crossings, not present in the simple 
Ising chain. The first of these is the transition between 
one F configuration of the Ising spins and the other: 
(tttt)---) <-> (iUI) ■ ■ ■)■ This transition takes place 
when £i(tt) = £i(44-)- This level crossing, induced by 
the magnetic field, can for example be encountered in 
frustrated Ising-Heisenberg chainsP^l A second kind of 
new ground state transition arises from the crossing of 
levels within a decorating unit. A ground state crossing 
can result in which the Ising spin configuration remains 
the same but the state corresponding to £i(s, s') crosses 
with the state corresponding to £ 2 (s,s'). More precisely 
this happens in the F configuration when £i(ft) = £ 2(tt) 
and in the AF configuration when either £i(fi) = £2(ti) 
or £i(|t) = £ 2(lt)i or both. The previous paragraphs 
have shown that we do not need to consider configura- 
tions other than F and AF for the ground state. 

Level crossings are usually connected with the presence 
of good quantum numbers. For the Ising-type crossings, 
the relevant conserved quantities are the n Ising spins 
{si}. The crossing of energy levels within the decorating 
unit should be associated with a conserved variable that 
is internal to that unit, much the same as in isolated 
molecules. In Section |III| we will encounter an example 
where both transitions — Ising type and internal type — 
occur in a magnetic field. 



D. Magnetization of powder samples 

In the following sections we will be comparing our the- 
oretical results with measurements performed on powder 
samples of the crystalline compounds. In this section we 
consider the powder averaging of magnetization for the 
example of the simple Ising chain. 

Let 9 and <j> be the polar angles of the magnetic field 
vector with respect to a molecular reference frame, then 
the free energy is a function of 8, (f>, and the strength of 
the field, B: f(6 1 <t> 1 B). The projection of the magneti- 



7 



zation on the field direction is 



e B -M(9,(f>,B) 



dB 



(21) 



Averaging over one hemisphere gives the powder magne- 
tization 

2tt Jo J dB 

Let us see how the powder averaging affects the magne- 
tization curve for a simple Ising chain. Take a spin-1/2 in- 
finite antiferromagnetic Ising chain with anisotropy axes 
parallel with each other and with the z axis, and uniax- 
ial g- factors (g x = g y = and g z = g). This could for 
example be realized by a chain of identical Dy 3+ units 



(see Section II B I . The Hamiltonian is given by Eq. ([lj , 
substituting 
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FIG. 3. Magnetization along the z direction versus field of 
the antiferromagnetic Ising chain, defined in Eq. ( 23 1 . The 
curve approaches a perfect step as T — ¥ OK. 



-JSiS 



flBgSiB COS 9, 



(23) 



where B cos 9 is the z component of the magnetic field, 
and J < 0. Defining j = J/k and b = /i^gB/k, the 
magnetization, which has only a nonzero z component, 
i^P 



M 7 . = 



fi B g sinh[6cos6 l /2T] 



sinh 2 [6cos(9/2T] +e~i/ T 



The projection on the field direction [Eq. (21)] is 
M z cos9. Plugging this in Eq. (22) and substituting 



u = cos 9 yields the powder magnetization of the Ising 
chain 



M 



sinh[6w/2T] 



sinh 2 [6w/2T] +e-J/ T 



-udu. 



(24) 



Because the Hamiltonian in Eq. (|23| does not depend on 



4>, this variable has been integrated out in Eq. (24). 

Figs. [3] and [4] show plots of magnetization versus mag- 
netic field, for coupling constant j = — 1 K. The step-like 
appearance of M z is associated with the ground state 
crossing that occurs at bcos9 — \j\. At that point, the 
antiferromagnetic ground state (or rather ground state 
Ising doublet) is replaced by the ferromagnetic state (all 
spins up). Consequently, the magnetization jumps from 
zero to the saturation value of 0.5, as seen in Fig. [3] 
The magnetization of a powder sample of the same Ising 
chain is shown in Fig. |4j Before the crossing point, M 
behaves qualitatively the same as M z . After the crossing 
point however, M is seen to reach only slowly its satu- 
ration value, which is half of the saturation value of M z , 
viz. 0.25. The limiting curve of M as T — > OK can be 



calculated exactly from Eq. ( 24 ) 



M 



T-s-0 




if b < \j\ 




FIG. 4. Magnetization of a powder sample of antiferromag- 
netic Ising chains. The low-temperature limiting curve is 
given in Eq. (25 1. 



(This is of course only valid for the antiferromagnetic 
case j < 0.) In Fig. |4j this limiting curve is very closely 
approximated by the curve at T = 0.01 K. Clearly, the 
sharp step of M z transforms in the powder to the concave 
form displayed by M. This is understood from the fact 
that, in a powder, for a given field b > there is always 
a fraction of molecules that is not magnetized (in the 
sense that they are in the antiferromagnetic ground state) 
because they are oriented so with respect to the field, that 
b cos 9 < j j I (see Fig. [3]) . The powder saturates only when 
every molecule is fully magnetized, and this happens only 
for b — > 00. Therefore, M (at OK) does not abruptly 
saturate at the crossover point, but increases slowly to 
saturation. 



E. 



Corrections for the contribution of excited 
Kramers doublets 



Heisenberg 
Ising 



In deriving the exchange Hamiltonian in Section |IIB| 
we assumed that only the lowest Kramers doublet on 
Dy 3+ took part. This is certainly a good approximation 
when the gap between the lowest and the second-lowest 
Kramers doublet is much larger then the strength of the 
exchange interaction. However, the excited Kramers dou- 
blets often have to be taken into account to a certain de- 
gree of approximation if a comparison with experimental 
data on susceptibility and magnetization is desired. The 
crystal field splitting of the 6 H 15 y 2 level is of the order of 
kT at room temperature. This gives rise to two effects: 
(i) a thermal population of excited Kramers doublets, 
and (ii) a modification of the lowest Kramers doublet as 
a function of the applied magnetic field by interaction 
with the excited doublets. 

Effect (i) is mainly visible in the temperature depen- 
dence of \T, where x is the powder magnetic suscepti- 
bility: for a single Dy 3+ center, %T increases monotoni- 
cally with increasing temperature, from the value of the 
ground doublet at K to the saturation value of 6 H 15 ^ 2 at 

higher temperatures. Effect (ii) gives rise to temperature- 
independent paramagnetism (TIP). It is visible at tem- 
peratures sufficiently low so that only the ground doublet 
is occupied. It contributes a linear increase of x^ with 
T and a linear increase of the magnetization M with the 
applied field B. 

In the simplest approximation, the contribution of the 
excited Kramers doublets to the magnetic properties of 
the chain is equal to the contribution they have to the 
properties of the single, isolated Dy 3+ ion in the same lig- 
and environment it has in the chain. Let xdic and Mdic 
denote susceptibility and magnetization derived from the 
decorated Ising chain model, and let XDy denote the sus- 
ceptibility of the Dy 3+ center and Mby-^ the magnetic 
moment induced by B in the ground doublet of the Dy 3+ 
center, then the corrected properties are (supposing one 
Dy 3+ ion per unit cell) 



(26a) 
(26b) 



XT = XDICT + [XByT - XByT\ T=Q/ , 

M = M mc • //„/>'. 



The last equation assumes that only the ground doublet 
of Dy 3+ is occupied. This is correct at the tempera- 
ture at which magnetization curves are usually recorded 
(e.g., 2K). If necessary, corrections due to other mag- 
netic ions can be added in the same way. XDy and 
fi'jjy can be obtained from multiconfigurational ab initio 
calculation d 7 * 8 * 22 *^, or, in oligonuclear complexes, exper- 
imentally by replacing certain magnetic ions with dia- 
magnetic ionsPi 1 



The equations ( 26 ) are evidently correct in the limit 
of vanishing exchange interactions. The assumption we 
make is that the exchange interactions are sufficiently 
small for these corrections to remain valid. A more ac- 




FIG. 5. Scheme of one unit cell (black) of the [DyCuMoCu] oc 
chain, showing type of exchange interactions and labeling of 
exchange constants. See Ref. [7J for the complete molecular 
structure. 



curate approach should take into consideration the fact 
that the excited Kramers doublets also participate in the 
exchange interaction with neighbors. This is however out 
of the scope of the decorated Ising model. 



III. [DyCuMoCuJoo POLYMER AS 
ISING-HEISENBERG CHAIN 

We now turn our attention to two actual examples of 
decorated Ising models based on Dy 3+ : a chain, treated 
in this section, and a four-ring treated in the follow- 
ing section. The problem is approached as follows: the 
Hamiltonian for the chain in a magnetic field is formu- 



lated, with the help of the considerations in Section II B 



Values of the g factors of the magnetic ions are taken di- 
rectly from ab initio calculations, reported elsewhere!^' 
This leaves the exchange coupling constants as parame- 
ters of the model, to be fitted by comparison with exper- 
imental magnetization and susceptibility data. (In Sec- 
tion |IV| the direction of the anisotropy axis contributes 
one extra parameter.) 

The [DyCuMoCu] chain was recently synthesized and 
details of its chemical composition and structure are 
given in Ref. [71 The crystal structure was found to consist 
of parallel linear chains each made of [DyCuMoCu] unit 
cells. Fig. [5] shows how the metal ions are connected by 
ligand bridges. Multiconfigurational CASSCF/RASSI- 
SO calculations have been performed on each of the four 
metal ions in their ligand environment, suitably discon- 
nected from the rest of the chain (details of the calcula- 
tions can be found in Ref. [7J and the accompanying Sup- 
plementary Information). Most important for us is that 
the Dy 3+ center was found to have a ground Kramers 
doublet, separated by 141 cm" 1 from the second dou- 
blet, and characterized by complete uniaxial anisotropy: 



Dy 



3+ 



g z = 19.6, g x = 0, g y = 0. 



(27) 



(Actually g x and g y were calculated about 0.03, which 
is small enough to be ignored.) The value of g z = 19.6 
shows that this doublet is only slightly perturbed from 
the \Mj = ±15/2) doublet of the 6 H 15/2 level, the lat- 
ter having g z — 20. Together with the fact that the 
energy gap to the second Kramers doublet is about ten 
times larger than the exchange interaction (as we will 
find later), these results indicate that the Dy 3+ ion will 



behave as an Ising spin, as described in Section II B As 



9 



a side-note we may add that the total splitting of the 
6 Hi5/2 level was calculated to be 560 cm -1 , which is in- 
deed of the order of room-temperature kT (see Section 
IlE l. 

Both Cu 2+ (d 9 ) and Mo 5+ (d 1 ) have a spin = 1/2, 
orbitally nondegenerate ground state, well separated (> 
15000 cm -1 ) from higher states. The two Cu 2+ ions in 
the unit cell reside in almost identical environments^ and 
have therefore virtually the same properties. The calcu- 
lated g factors are tetragonal: g\\ = 2.33, g± = 2.07 for 
Cu 2+ and , 9 || = 2.00, g ± = 1.95 for Mo 5+ . To avoid 
unnecessary complications we will regard these ions as 
isotropic spins with root-mean-square g factors 



Cu 2 
Mo 5 



ffcu = 2.16, 
9mo = 1-97. 



This approximation will not have important conse- 
quences for the magnetic properties, which are largely 
dominated by the high Dy 3+ moment anyway. 

We introduce exchange interaction between metal 
ions directly connected by ligand bridges. Dy 3+ in- 
teracts with its three neighbors via the Ising Hamilto- 
nian Eq. (14 1. Interaction between the isotropic spins is 
given by the Heisenberg Hamiltonian — JjjS, • Sj. Fig. [5] 
shows the exchange configuration, with single bonds rep- 
resenting Ising interaction and double bonds representing 
Heisenberg interaction. Note that we have approximated 
the Dy-Cui and Dy-Cu2 coupling strengths to be equal 
(Ji), following the approximate local symmetry of the 
Dy-Cu pairsP 

It is now possible to see that the [DyCuMoCuJoo 
polymeric chain is indeed an experimental realization of 
an Ising-Heisenberg chain where the [CuMoCu] trimeric 
Heisenberg units decorate the Dy-Dy bonds and the 
Dy 3+ Ising spins separate the [CuMoCu] Heisenberg 
trimers from each other. The total Hamiltonian in a 
magnetic field B is then given by Eq. and 

^i( s ii s i+l) = — ^l s i<SjCui — JzSiS- Mo — J3S1CU1 ' SiMo 
— MBffDySi-Bz — MB(3CuSjCui + <7MoSjMo 

+ 3cuS 4 cu 2 ) • B, 

(28) 

where sf is shorthand for sf Dy and guy is the g z factor of 

Dy 3+ [Eq. p7[)]. The z axis is the anisotropy axis of the 
Dy 3+ center. We have not specified its direction with 
respect to the chain axis but this is not important here 
because there are no other axes in the problem (the Dy 3+ 
anisotropy axes are parallel by translational symmetry 
and we have assumed Cu 2+ and Mo 5+ isotropic). All 



the spins in Eq. (28 1 are spins of 1/2. 

The Hamiltonian exhibits some symmetry. It is ro- 
tationally invariant around z, if B is rotated simulta- 
neously. We may therefore restrict B to lie in a plane 
through z, say the xz plane. This simplifies calculation 
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FIG. 6. Powder magnetization of [DyCuMoCu], 



of the powder magnetization Eq. (22): one has to inte- 



grate only over 8. When B is directed along the z axis, 
the z component of the total spin in each decorating unit 
is conserved: 



S~ — S, 



iCui 



-s 



iMo~ 



-s 



iCu 2 ' 



[S?,H(B\\z)}=0. (29) 



We also note that in this case the Zeeman Hamiltonian 
commutes almost with H. It would commute exactly 
when gcu — 5Mo, for then the last term in Eq. (28) re- 
duces to —fiB9CuSfB z . 

We let the length of the chain go to infinity: n — > oo. 
To solve for the thermodynamic properties we are only 
require d to fi nd the eigenvalues £fc(s, s') of Eq. ( p8| (See 
Section II A I, with k = 1 . . .8, corresponding to the 2 3 



possible states of the [CuMoCu] spin unit. This is done 
by 4 numerical 8x8 matrix diagonalizations, one for each 
(s, s') pair. 

We can now compare the theory with experiment. 
Powder magnetization (at 2 K) and susceptibility data 
have been recorded.^ We recall that we have to cor- 
rect the theoretical curves before comparing with exper- 
iment according to Eq. (26). The corrections are pro- 
vided by the ab initio calculations.^ /Lt^ turns out to be 



0.03 /ib/T; the accompanying correction in Eq. (26b) is 



never more than 2.5% of M. We ignore this correction. 



We do however correct x as in Eq. (26a 



Closest agreement with experiment was found for the 
following values of the exchange constants (plots in 
Figs. § and |7|): 



Ji = 15.3 cm , J 2 
Ja = 11.8cm" 1 . 



5.0 cm" 1 , J 3 



-8.3cm" 1 , 

(30) 



These were obtained by a least-squares fit of x followed 
by a small manual adjustment to improve the fit of M 
while not distorting that of x appreciably. (The least- 
squares fit of x gave J\ = 15.7cm" 1 , J 2 — —8.3cm" 1 , 
J3 = —6.3 cm" 1 , J4 — 11.8 cm" 1 .) There are, as far as 
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FIG. 7. Powder magnetic susceptibility of [DyCuMoCu]oo- 
The theoretical curve contains the correction for the contri- 
bution of excited Kramers doublets. 
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FIG. 8. Powder magnetic susceptibility of [DyCuMoCu]oo. 
The correction refers to the contribution of excited Kramers 



doublets, Eq. (26a 



we know, no data in the literature with which to compare 
the values in Eq. ( 30 ). However, an experimental study is 
available of a (Dy , Cu 2+ ) dinuclear complex in which 
the bridging ligand is the same as in this chainPS The 
authors found a ferromagnetic interaction. A superficial 
analysis of the susceptibility curve in that paper, using 
the Ising Hamiltonian we use in this paper, yields J\ — 
15 ± 5 cm" 1 . The other values in (30 1 arc difficult to 
assess. For a discussion of these values and their relation 
with the molecular structure as well as some evidence 
from DFT calculations, we refer to Ref. Certainly, no 
confidence should be attached to the numbers in decimal 



places in (301 



The effect of the excited Kramers doublets of Dy is 
most clearly seen in the curve (Fig. [8}. The curve 
shows a steady increase above 50 K which is not pre- 
dicted by our decorated Ising model, but is indeed due 
to the thermal population of the Kramers doublets that 
originate from the 6 H 15 ^ 2 level. We can obtain the ex- 



pected high-temperature limit of \T by considering the 
metal ions as independent spins. The susceptibility com- 
ponents Xaa of an angular momentum multiplet J with 
principal g-factors g a (a = x, y, z) are given by^ 



Xaa 



3fcT 9aJ(J 



Summing over Dy 3+ ( b H 15/2 , g 15/2 
and Cu 2 (all are isotropic) gives 



1). (31) 
= 4/3), Cu 1; Mo, 



XT = - Xaa 



T 



Na^ b 
3k 



015/2' 



15 17 

2~T 



(9 2 cu 



2 

9mo 



9cu) 



1 3 
22 



= 15.4- 



mol 



A similar calculation, only including the lowest Kramers 
doublet of Dy 3+ , with g- factors as in Eq. (271, gives 
13.2 cm 3 Kmol -1 . The correction supplied by the ab ini- 
tio calculations to account for this difference, is seen to 
cover nicely the high-temperature part of the experimen- 
tal curve. 

One notices that xT shows a slight depression around 
40 K which is not entirely reproduced by the theory. This 
might indicate a failure of the simple approximation we 
used to include the excited Kramers doublets. Eq. (26a) 
is certainly correct at very high temperatures, when the 
exchange interactions are irrelevant, and at very low tem- 
peratures, when the excited Kramers doublets are not 
occupied. If these two regions do not overlap, however, 
there is a temperature window between, in which excited 
doublets start to get occupied while exchange interaction 
is not quite negligible yet. In that case, the exchange in- 
teraction of the occupied excited doublet (s) with other 
ions should be taken into account. Such an interaction 
of antifcrromagnetic type could possibly depress xT as 
observed. 

We shall now describe some features of the spectrum of 
the chain, paying attention to the properties described in 
Section [II C| Consider the chain without magnetic field. 
The exchange parameters in Eq. ( 30 1 predict a ground 



state that has an AF Ising spin configuration. This is in 
accordance with the susceptibility measurements, which 
show that xT -> as T -> 0, requiring a nonmagnetic 
ground state (Fig. 8]). The ground state is indeed non- 
magnetic because fei (T-i)} is the time-reversed state of 
| £l (;t)>. Let M s 
{- 



denote an eigenvalue of S z [Eq. ( 29 )] : 
4/2,1/2,3/2}. Sf is conserved so M s 



M s G {-3/2, 
may be used to label the eigenstates \sk(s,s')) (we may 
leave out the index i because all units of the chain are 
identical). For the ground state, we find Ms = —1/2 in 
|£i(ti)) and M s = 1/2 in |e a (4-t)>- 

Since the ground state is AF, we might expect that in 
a magnetic field a crossover will occur to an F ground 
state. This is indeed what happens. The convex increase 
of M in Fig. [H] points to a flip of the Dy 3+ spins to a 
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FIG. 9. Theoretical magnetization of [DyCuMoCu]oo- The 
powder magnetization (See also Fig. [fi]| is compared with the 
projections of M on the field direction [Eq. (|2 1 p] , for three 
different directions of the field; 8 is the angle between B and 
the z axis. 
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FIG. 10. Theoretical magnetization of [DyCuMoCu]oo. Same 
as Fig. [9] but at lower temperature and to higher field. 



parallel configuration. This is inferred from the value of 
the magnetization, which approaches 6 /iB at 5 T. The 
[CuMoCu] unit alone can only contribute a maximum of 
(2.16+1. 97+2. 16)/2 = 3.15 /j, b . The strong increase must 
come from the contribution of the large Dy 3+ moments. 

The behavior of magnetization along certain directions 
of applied field is shown in Fig. [9j The AF — > F transition 
is most clearly seen when the field is applied along z 
(9 = 0); the transition occurs below IT. After IT, M z 
reaches an approximately constant plateau at ~ H/iB- 
The saturation value of magnetization in direction < 
9 < tt/2 is (19.6 cos 6 + 2.16 + 1.97 + 2.16)/2. This gives 
12.9 /iB for 9 = 0, which shows that M z has not quite 
reached its maximum at 5T yet. 

The positions of level crossings become more sharply 
defined on lowering the temperature (Fig. 10). Here we 
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FIG. 11. Eigenvalues of h(s, s') [Eq. (28l] in a magnetic field 
parallel with the z axis (8 = 0). Circles indicate ground 
state level crossings. The ground state of the chain is AF 
in zero field (left), switches to F at 0.64 T, and undergoes 
an internal level crossing at 6.3 T, marked by a change of 
the internal quantum number Ms from 1/2 to 3/2. Both 
crossings can be seen in the 9 = magnetization curve in 
Fig. |10| Note that the energy curves appear as straight lines, 
although, with the exception of Ms = ±3/2, they are not 
exactly straight, because the Zeeman Hamiltonian does not 
completely commute with the total Hamiltonian. All £fc(tt 
) decrease with increasing field strength because the large 
magnetic moment of Dy 3+ dominates the smaller magnetic 
moments of the decorating unit. 



also see that M z undergoes a second transition at 6.3 T, 
after which it reaches saturation. This transition is con- 
nected with a level crossing in the [CuMoCu] unit (see 
Section II C I from Ms = 1/2 to A/5 = 3/2, as opposed to 



the first transition, at 0.64 T, which is of the Ising type, 
described by Eq. ( 20 ) . The latter is the analogue of the 
transition in the AF simple Ising chain (Fig. [3]), while the 
"internal" transition has no such analogue but is unique 
to the decorated Ising chain. The relevant energy level 
diagram is shown in Fig. [IT] Note that, for fields not 
parallel to z (for example, 9 = 7r/4 in Fig. 10 1, Ms is 



not a quantum number and the internal level crossing 
turns into an avoided crossing. This does not apply for 
the Ising level crossing because the Ising spins are always 
conserved. Only when the field is applied perpendicular 
to z (6 = tt/2 in Fig. 10 1 does the AF — > F transition 
not occur because the Dy 3+ spins do not interact with 
perpendicular fields. 
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FIG. 12. Theoretical susceptibility of [DyCuMoCu]^, with- 
out correction for contribution of excited Kramers doublets. 
The powder yT (see also Fig. [8| is compared with the Carte- 
sian components of yT . z is the direction of the anisotropy 
axis of Dy 3+ , x is any direction perpendicular to z. X ~ 
(Xz* +2xxx)/3. 



The low-temperature limit of the powder magnetiza- 
tion in Fig. [lO] may be compared with that of the sim- 
ple Ising chain in Fig. [4] The resemblance is clear; the 
decorated chain is different in the small linear increase of 
M before the transition, and the more linear approach to 
saturation, which lies at (19. 6/2 + 2. 16 + 1.97 + 2. 16)/2 = 
8.0 /ib- Both are due to TIP interaction in the [CuMoCu] 
unit, the effect of which is most clearly seen in the M x 
curve in Fig. [TU] 

To conclude this section we remark that the mentioned 
similarity with the magnetization of the simple Ising 
chain is a consequence of the very high magnetic moment 
of the Dy 3+ spins in comparison with the [CuMoCu] unit. 
The dominance of Dy 3+ is most dramatically shown 
in the components of \T (Fig. 12). An application of 
Eq. (|31[ ) shows that the high-temperature limit of XxxT is 

^a(2.16 2 + 1.97 2 + 2.16 2 )i| = 1.24 cm 3 Kinol" 1 , while 

that of XzzT is ^a(19.6 2 + 2.16 2 + 1.97 2 + 2.16 2 ) 
37.2cm 3 Kmor 1 .' 
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IV. Dy 4 Cr 4 COMPLEX AS DECORATED ISING 
RING 

As a second example we describe in this section the ap- 
plication of the decorated Ising model to a ring-shaped 
Dy 4 Cr 4 moleculeP Dy 4 Cr 4 consists of alternating Dy 3+ 
and Cr ions forming a closed ring. The four Dy 3+ ions 
lie in a plane. The Cr 3+ ions are positioned alternat- 
ingly above and below this plane, "decorating" the Dy- 
Dy bonds. The molecule has Z?2d symmetry, the Dy 3+ 
ions lying on C2 axes and the Cr 3+ ions lying on the 
mirror planes. We choose a molecular reference frame 
XYZ so that Z coincides with the S4 axis and X and Y 
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FIG. 13. Schematic representation of the Dy 4 Cr 4 molecule 
indicating numbering of atoms and exchange coupling con- 
stants. The boxes show the orientation of the local anisotropy 
axes on Dy sites, when viewed from the poles of the X and 
Y axes. The Z axis points out of the center of the scheme. 



Ab initio calculations have been performed in the same 
way as for the [DyCuMoCu] chainJS From these, we take 
again the g-factors of the ground doublet of Dy 3+ and 
of the isotropic ground state spin multiplet of Cr 3+ (d 3 , 
5 = 3/2): 



3+ 



Dy 
Cr 3 + 



g z = 19.66, g x = 0, 
9Cr = 1-97. 



9y 



0. 



(32) 



The Dy 3+ Kramers doublet is again very close to the 
\Mj = ±15/2) state, permitting the use of the Ising 
model. However, the same calculation predicted the sec- 
ond Kramers doublet at 30 cm -1 , not very high com- 
pared with exchange interaction, which we found in the 
previous section w 10cm _1 . This should be seen as a 
warning that our treatment of the excited Kramers dou- 
blets as "innocent" may not be entirely correct here, 
and thus may lead to discrepancies with experiment. 
In this respect we must also note that, for such small 
excitation energies, the results of the ab initio calcu- 
lations are not always conclusive on the nature of the 
ground state Kramers doublet. In the present case, for 
instance, another set of calculations produced a ground 
state Kramers doublet on Dy 3+ that is not uniaxial as in 
Eq. (32 1, having relatively large transversal g-factorsP 



The decorated Ising model would be unusable in this 
case. We find however that, assuming the axial g- factors 
in Eq. (32 1, Dy 4 Cr 4 is an interesting example of a dec- 



coincide with the two C2 axes of Z?2d (Fig. 13 ) 



orated Ising ring, for which qualitative agreement with 
experimental magnetic properties can be obtained. 
Direct ligand bridges connect each Dy 3+ with two 



neighboring Cr 3+ ions and two neighboring Dy J+ ions 



3+ 
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Exchange interaction between these pairs is introduced 
[Eqs ( fl3| and (14)]. The Hamiltonian is then given by 
Eq. MTh = h!+h 2 + l 



h-2 + /13 + h,4, and 

j / Zi ryZi 1 2 i + l r< z i + l \ T Z% 2 

-Ji(V<V + Si ) - ■hs i z s i 

■ UBig-DySi' B Zi + gCrSi ' B), 



+ 1 
hi 



(33) 



where s Zi denotes the Ising spin- 1/2 variable on Dy^ and 
S Zi denotes the projection of the spin of Cr^ on the mag- 
netic anisotropy axis of Dy^ (for numbering, see Fig 13). 
Similarly, B Zi is the projection of the magnetic field on 
the anisotropy axis of Dy.;. gD y is the g z factor of Dy 3+ 
[Eq. pi]. 



An interesting difference with the [DyCuMoCu] chain 
is that here, in Dy 4 Cr 4 , the four anisotropy axes are 
not, in general, parallel, a result of point symmetry in- 
stead of translational symmetry. The orientation of the 
local anisotropy axis on Dy 3+ , being one of the <?-tensor 
main axes, is restricted by the local C 2 symmetry to be ei- 
ther parallel with, or orthogonal to the local C 2 axis. The 
first possibility can be excluded on the basis of the ex- 
periment; with the Zi pointing radially outwards at each 
Dyj , the ground state of the whole molecule is necessarily 
nonmagnetic, because the local moments add up to zero, 
independent of whether the ground state is F or AF with 
respect to the Ising spins. The experimental suscepti- 
bility measurement however indicates a magnetic ground 
state (nonzero intercept on the vertical axis in Fig. 14). 
We must therefore choose the second case and let the 
anisotropy axis on each Dy be orthogonal to the local 
C2 axis and make an angle of a with the molecular Z- 
axis (See Fig. 13). By applying the symmetry elements 
of D 2 d to one of these anisotropy axes, one obtains the 



other three. When a = the four axes are parallel and 
point in the same direction as Z. We note that the ab 
initio calculations yielded a = 37°. We will need some 
flexibility in our model however, so we leave a as a pa- 
rameter that will be determined from comparison with 
experiment. 

In terms of the molecular coordinate system, the pro- 
jections on the local anisotropy axes are a function of 
a: 



Sf = cos(a) Sf + sin(a) Sf 
Sf = cos(a) Sf + sm(a)Sf 
S Z3 = cos(a) Sf - sm{a)Sj 
= cos(a)5f - sin(a)5f . 



(34) 



The same relations hold for the magnetic field, after re- 
placing Si by B. 

The fact that only exchange interactions of Ising type 
appear in Eq. ( 33 ) makes it possible to find analytical 



solutions of the eigenvalues and the partition function. 
From Eqs. (33) and (34) we see that the part of hi that 



involves S% is a projection of Si on the vector 

- Ji(s 4 z, + s i+ iz i+ i) - gcrMBB, (35) 
where Zj is the unit vector along the anisotropy axis of 
Dyj (the superscripts Zi on Si are left out from now on). 
The vector ( 35 1 defines the quantization axis of Sj , which 



depends on the states on the neighboring Dy + sites 
(sj, Si+i). The stronger the coupling (Ji) with Dy, the 
stronger will be the deviation of the quantization axis 
from the direction of B. The eigenvalues of hi are then 

e iM s(si,Si+i) = hM s - J 2 8is i+ i - iXB9v y SiB Zi , (36) 

where Ms — —S, . . . , S, and 



I 2 

T (1 



AsiSi+i cos 2 a) + 2J 1 [i B gcr(s i B z - + s i+1 B z ^) + ^g 2 Cr B 2 



(37) 



is the length of the vector in Eq. ( 35 1 . Some remarks 
should be made on the solutions. Eqs. (37) and (34) (re- 



place Si by B) show that the spectrum in Eq. ( 36 ) is not 
the same for every unit i, as it was in the [DyCuMoCu] 
chain, unless B is applied along the Z axis. This means 
that also the transfer matrices Tj will be different and 
that we have to use Eq. ([7| instead of Eq. (|8| for the 
partition function. A second remark concerns the quan- 
tum number Ms- The lowest energy in Eq. (36) is al- 
ways given by Ms = —S, but note that the axis to which 
this quantization refers is not invariant; in particular, it 
changes with strength and direction of applied field, so 
that Ms does not represent a real conserved quantity that 
could be responsible for level crossings of the "internal" 
type. Such crossings do not occur in Dy 4 Cr 4 . 



We conclude the solution by finding the partition func- 
tion Z, Substituting Eq. Q in Eq. Q we find 

s 



E 



M s =-S 

sinh[/3& i (2S + l)/2] 



sinh[/3&i/2] 
x exp[(3(J 2 SiSi +1 + HBgr> y SiB Zi )} 

With Tj as defined in Eq. we obtain the partition 
function 

Z = Tr(TiT 2 T 3 T 4 ). 

Let us now compare the theoretical results with experi- 
ment. A great amount of information on the values of the 
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14 1. The nonzero inter- 
indicates a magnetic 



parameters a, J\ and J2 can be obtained by inspection 
of the powder yT curve (Fig. 
cept xT| T ^ = 34.6 cm 3 K mo 

ground stated Now from the general theory we know 
that the ground state is either F (tttt) or AF (tlti-) 
with respect to the Dy 3+ spinspSl AF is nonmagnetic so 
we decide that the ground state must be F. 

Incidentally, we can precisely delineate the regions in 
parameter space where the ground state is F or AF: 



F (fttt) : J2 > -V2S (Vl + cos 2 a-sma) (J^, 

(38a) 

AF(tlU) : J2 < -V2S Ufl + cos 2 a- sin a) \J X \. 

(38b) 

A second piece of information comes from the increase 
of yT with increasing temperature. This is partly but 
not completely due to the occupation of excited Kramers 
doublets, as one can show by subtracting the contribu- 
tion of the latter, obtained from the ab initio calculations 
(not shown here). There must still be an antiferromag- 
netic interaction to explain the increase. Since the Dy 3+ 
are already known to be ferromagnetically aligned, the 
only possibility is that the Cr 3+ spins couple antiferro- 
magnetically with Dy 3+ , or J\ < 0. 

With this information, we can determine the angle 
a. At OK, \T is determined by the magnetic moment 
in the ground state onlyP32D l n the F (tttt) state, 
XxxT\ T -+o = XyyT\ t ^o = by symmetry and 



XzzT\t^o 



4 4 

# D y COS a ^ s i+9Cr ^2 
i=l i=l 



Sf\F 



(39) 

so X T\t^o = \xzzT\t^o. With the help of Eqs. §4|- 
(37) and the fact that, in the ground state, Ms = —S in 
Eq. |36| , we can evaluate Eq. ( 39 1 to find 



XT\ 



T-i-0 



+ sgn(Ji)2 V / 2fl'Cr5 



+ cos 2 

This is a strictly decreasing function of a that can be 
used to derive a from the experimental value xT\t^q = 
34.6cm 3 Kmol _1 , and the knowledge that sgn(Ji) = —1. 
This gives a = 49°. 

We derive values for Ji, J2, and a by a least-squares 
fitting of x- As before, a correction for the contribution 
of excited Dy 3+ Kramers doublets is provided by the ab 
initio calculations and applied following Eq. ( 26 1 . The 
fitting yields 



Ji = -8.9 cm' 



Jo = 5.2 cm 



a = 49.2°. 



The comparison with experiment is shown in Figs. [14] 
and |15| Note that the magnetic properties are reported 
per mole (x and %T) or per molecule (M) of Dy 4 Cr 4 and 
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FIG. 14. Powder magnetic susceptibility of Dy 4 Cr 4 . The 
correction on the theory refers to Eq. ( 26 1 . 
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FIG. 15. Powder magnetization of Dy 4 Cr 4 . A linear correc- 
tions of 0.5/xb/T has been added to the theoretical curve, 
according to Eq. (26b I. 



not per DyCr unit. Note also that sgn(Ji) = —1, that 
J2 satisfies Eq. (38a), and that a agrees with the value 
derived above. 



The agreement of magnetization curves (Fig. 15 1 is not 
as good as it was for the [DyCuMoCuJoo chain, although 
the qualitative properties seem to correspond. In par- 
ticular, we mention the strong linear increase of M at 
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higher fields (~ 2.6/is/T), which is due to the gradual 
orientation of the Cr 3+ spins to the magnetic field [see 
discussion connected with Eq. (35)], and, to a smaller ex- 
tent, also to the correction of 0.5/ie/T, a non-negligible 
linear contribution to magnetization, which is due to the 
low-lying excited Kramers doublets. 

As was mentioned before, the discrepancies are not un- 
expected given a low-lying first excited Kramers doublet 
of Dy 3+ , which could undermine the assumptions under- 
lying the decorated Ising model. Note also that we could 
not take the ab initio value of 37° for a. Leaving a as a 
parameter can be seen as a partial compensation for the 
inaccuracies of the model and the ab initio results. 



V. CONCLUSION 

We have shown that the decorated Ising model is 
a valid model for the magnetic properties of certain 
lanthanide-containing magnetic compounds, if the crys- 
tal field spectrum of the lanthanide ion satisfies certain 
properties. The most important of these is the require- 
ment of a ground state Kramers doublet with completely 
uniaxial magnetic anisotropy (this statement is simpli- 
fied, see Section II B for the correct details). It is a re- 



markable fact that precisely this property has been estab- 



lished by multiconfigurational ab initio calculations on 
several Dy 3+ centers that are part of polynuclear molec- 
ular magnets. Perhaps the best known example is the 
Dy 3 triangle, where the Ising properties of Dy 3+ were 
used to explain the nature of the ground statePZI 

We have focused on Dy 3+ as lanthanide ion because 
this is a much-used lanthanide in current synthetic re- 
search in molecular magnetism (witness both compounds 
in this paper) and because computational results show- 
ing that it meets the requirements for an Ising spin are 
available. However, there is no reason to assume that the 
findings are unique to Dy 3+ . We expect that other lan- 
thanides with high momentum (e.g., Er 3+ ) will exhibit 
the same uniaxial anisotropy in certain ligand environ- 
ments and that examples of decorated Ising chains based 
on lanthanides other than Dy 3+ will be found in the fu- 
ture. 
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